Non Poisson intermittent events in price formation 
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The formation of price in a financial market is modelled as a chain of Ising spin with three 
fundamental figures of trading. We investigate the time behaviour of the model, and we compare 
the results with the real EURO/USD change rate. By using the test of local Poisson hypothesis, we 
show that this minimal model leads to clustering and "declustering" in the volatility signal, typical 
of the real market data. 
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Despite the obvious interest to adequately describe 
the stochastic behaviour of price formation in a financial 
market di] , stochastic modelling misses some interesting, 
and perhaps crucial characteristics of the real market. In 
particular, it is well known that the market cannot be 
adequately described by time series of independent and 
identically distributed realizations of a random variable. 
In fact time series in the financial framework depend on 
a large number of strongly interacting systems. Then, 
the problem should be investigated in the framework of 
complex systems. In this case long-range correlations, 
not obvious a priori and very difficult to be adequately 
modelled, are present in the system. In this situation 
every step towards the goal of analysing and describing 
correlations is welcomed. The Ising spin system 2] is one 
of the most frequently used models in statistical mechan- 
ics. Being made simply by binary variables, the model is 
able to reproduce different complex phenomena in differ- 
ent areas of science like economy 0, 0, , biology @ and 
sociology , were many interactions are required among 
various discrete elements. Some times ago a model of 
this kind has been used to describe the decision making 
mechanism of a closed community 8j, where, in contrast 
to the usual majority rules |9(, the influence was spread- 
ing outwards from the center. 

In this paper we use the Ising spin as a minimal model 
of financial market, in order to investigate and repro- 
duce the non-Poisson characteristics of isolated events in 
price formation. Besides the statndard statistical tools, 
we use a local analysis of time intervals between volatil- 
ity events, which is not affected by the non-stationarity 
of the signal. 

The dataset we investigate is the exchange rate 
EURO/USD collected every minute during two months 
in the period from May 25, 2004, up to July 31, 2004. 
The exchange rate e(t) is represented by a discrete time 
series where t = tk = kAt (At is the sampling rate 
and k = 0,1,..., TV). The main quantities of interest 
for traders within the market is the return of price r(t), 
defined as the change of the logarithmic price over a time 
interval r, namely G(t) = In e(t + r) — lne(f). A measure 
of the magnitude of fluctuations of prices used by traders 
is the volatility, defined as the absolute values of returns 



averaged over a window of time extension T = nAt 

. t+n-l 

VT(t) = - £ |r(f)| 



(1) 



where n is an integer representing the moving average 
window size. A further quantity is the autocorrelation 
function of returns (or of volatility) defined as 
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where G n = G(t n ). The main characteristics of the time 
serie are summarized in Figure where we show the time 
evolution of various quantities, namely e(t), G(t) and 
Vr(t). We note the typical non-stationarity, the pres- 
ence of bursty and explosive events in the return dataset 
and of isolated volatility events. 

We also report the autocorrelation function for both the 
returns and the volatility (see Figure |2J). While the for- 
mer displays a very fast decay (m T < 10 min), the latter 
displays a much longer persistent time lag. 
Finally in Figure [3] we report the standardized probabil- 
ity density function (PDF) of returns, which shows very 
clearly typical "fat" tails. 

The above characteristics can be adequately modelled 
by the Ising spin chain. The community of interacting 
agents which partecipate to the market, is modelled as 
a chain of one-dimensional Ising spins Sj (i = 1, . . . , N). 
Each agent, at a time t, has an attitude to place buy 
orders, described by an up-spin Sj(i) = 1, or to place 
sell orders which is described by a down-spin Si(t) = — 1 . 
The price of the asset at a given time t is proportional to 
the normalized difference between demand and supply 



v(t) 



1 N 



(3) 



which in the classical Ising model is the magnetization. 
Here \y(t)\ < 1 represents a probability. The evolution 
of the Ising chain is controlled by three simple rules that 
describe three figures partecipating to the market as in 
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FIG. 1: We report the time evolution of the exchange rate 
EURO/USD e(t) (upper panel), the return G(t) with At = 
lmin (middle panel) and the volatility Vr(t) for T = 8min 
(lower panel). 
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FIG. 2: The solid line is the autocorrelation function for re- 
turns of the exchange rate EURO/USD displayed in lin-log 
scale. The dashed one represents the autocorrelation func- 
tion for the volatility of the same signal in the same scale. 

Ref. @: 

1) in general a lot of market agents follows the trend. 
This is described by the rule that both agents in two next 
neighbouring sites (i — 1, i + 2) follow the trend given by 
agents placed in (i, i + 1), namely when s^Si+i = 1 then 

Si-l = Sj+2 = Si. 

2) Some agents prefer to act on the market independently 
on the neighbour agents. Then in absence of a precise 
"tendency", sites (i — l,i + 2) are chosen at random, 
namely when SiSi + i = —1 then both Sj_i and Si + 2 are 
randomly chosen in an independent way. 

3) A third figure is partecipating to the market, that is an 
agent with some knowledge external to the market. This 
is known in literature as "fundamentalist" 0, . This 
figure has an exact knowledge of what happens in the 
market so that he knows the difference between demand 
and supply. If supply is greater than demand he places a 



FIG. 3: Standardized PDF for returns G(t) of the 
EURO/USD exchange rate. The solid line is the stretched 
exponential fit with with c = 1.0 ± 0.1. 

buy order, while he places a sell order when the demand 
is greater than supply. 

Realizations of the Ising spin model have been analysed 
starting from a random sequence of N = 1000 spins at 
t = 0. Each interaction between two agents is interpreted 
as the tick-by-tick interaction between traders, while the 
price fixed at a given time t, defined as p(t) = exp [j/(£)], 
is calculated after a single evolution of the whole chain. 
Of course time is arbitrary in the Ising model. Actu- 
ally, we have varied the value of the sampling rate in 
the numerical dataset and we have compared it with the 
real dataset. On the basis of the results obtained fron 
the simulation, we can suppose that each update of the 
whole chain represents the price every minute of the real 
market. Then, in the following, p(t) is proportional to 
e(t). 

As in Ref. Q , we compare the real data of EURO /USD 
exchange rate with the numerical results of the simula- 
tion. The statistical analysis of p(t) shows the same fea- 
tures of real exchange e(t), namely the non-stationarity 
of the signal (see Figure , the behaviour of the corre- 
lation functions (see Figure 0) and, finally, the fat tailed 
PDFs of returns, as shown in Figure HO 

The statistical analysis carried on until now shows that 
this minimal model of financial market can capture many 
features similar, in a statistical sense, to the real market 
ones. 

To further characterize the temporal properties of these 
time series (exchange rate and stock price) and investi- 
gate on the characateristics of isolated events in pricing 
formation, we analyze the distribution of waiting times 
(WT) between these events. 

The statistical distribution of WT between discrete or 
isolated "events" has been shown to be a powerful tool 
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FIG. 4: We report the time evolution of the stock price p(t) 10 1 1 ' 

(upper panel), the return G(t) with At = 1 (middle panel) -4 -2 2 4 

and the volatility Vr(t) for T — 10 (lower panel). 
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FIG. 5: The solid line is the autocorrelation function for re- 
turns of the stock price displayed in lin-log scale. The dashed 
one represents the autocorrelation function for volatility of 
the stock price in same scale. 

to investigate and characterize temporal point processes 
like solar flares [13, earthquakes [Tl.j, dissipative events 
in turbulent flows |l2j , disruptions in laboratory plasmas 
[l3| and so on . We consider the volatility of a finan- 
cial market as a temporal point process with a certain 
number of discrete "events" , as the spikes in Figures Q] 
andU 

To quantify these events, we need to fix a threshold, 
above which we can select a spike as an "event". We 
adopted a method according which all events that are 
above a certain threshold are cut. Then, we calculate the 
flatness of the remaining signal. The value of the thresh- 
old is varied, until the flatness of the remaining signal is 
equal to the Gaussian value 3 ^5|. At this point we cal- 
culate the waiting times from the difference between the 
end time and start time of two next events (quiet times). 
The statistical analisys of such variables can be used 
to detect the presence of long-range correlations in the 



FIG. 6: Standardized PDF for returns G(t) of the stock price. 
The solid line is the stretched exponential fit with with c = 
0.85 ±0.01. 

dataset. For example, the study of the autocorrelation 
function of quiet times of long-term persistent records 
can provide informations about the presence or not of 
clustering in the extreme events [l^ EH- D ue to the 
non-stationarity and the poor statistics of our dataset 
(~ 300 events) this kind of analisys can not be per- 
formed. However, a local Poisson hypothesis test can 
be used to check for presence of clustering in short and 
non-stationary records |171 Il8| . 

Let us consider the time sequence of extreme events 
selected in the volatility data according the rules ex- 
plained above. Let us suppose that each event occurs 
at the discrete time tf, each ith event can be associ- 
ated to the two nearest and next nearest neighboring 
events, which occur at the times £»— i, U+i, U-i and 
fi+2, respectively. Let SU = min{tj+i — UjU — U-i} 
and let #Tj be either = ti+2 — U+i (if SU — U+i — U) or 
Sri = U-i — U-2 (if Sti = ti — ti-i). If a local Poisson 
hypothesis holds, both SU and 5^ are indipendently dis- 
tributed according two Poissonian probability densities: 
P(Sti) — 2Xiexp(—2XiSti) and P(STi) = Xi exp(— XiSri) . 
In this case, the stochastic variable H defined as 



SU + \Sn 

is uniformely distributed in [0; 1] and it has a cumulative 
distribution P(H > h) = 1 — h independent on A. In a 
process, where SriS are systematically smaller than 2<5tjS, 
much more events are found respect to the Poissonian 
case, that is clusters are present and the average value 
of H is greater than 1/2. Conversely, a value of (H) 
less than 1/2 means that voids are present in the signal. 
The values of H have been calculated from the volatil- 
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FIG. 7: In the panel (a) probability P(H) of exchange rate 
EURO/USD volatility (solid line) and of stock price volatil- 
ity (dashed line) as a function of H is shown. A uniformely 
distribution P(H) in [0; 1] is also shown as a reference (dot- 
ted line). In the panel (b) cumulative probability P(H > h) 
of exchange rate EURO/USD volatility (solid line) and of 
stock price volatility (dashed line) as a function of h is shown. 
P(H > h) = 1 — h is also shown as a reference (dotted line). 



for the events obtained from both dataset. Indeed, the 
Kolmogorov-Smirnov test gives a probability of 1.7% for 
Poissonian distributions of WT of both dataset. The av- 
erage value (H) = 0.47 ± 0.3 and the shape of P{H) 
in Figure Ufa) indicate that voids and clusters are both 
present in a same fashion in both signals. 
Several observational evidences of clustering in the mar- 
ket volatilities can be found in literature and numerical 
models El illll. They refer to the clustering volatility 
phenomenon as a clustering of time periods when volatil- 
ity fluctuates very strongly and its autocorrelation func- 
tions show very long persistent time lags. Several ap- 
proaches have been guessed to describe this phenomenon, 
as in Ref. [22 • Here, by using the test of local Poisson 
hypothesis for the first time, we can provide a quanti- 
tative measure of the phenomenon of volatility cluster- 
ing in financial time series and we show that voids are 
also present. In summary, we have shown that a mini- 
mal spin model leads to a fat-tailed distribution function 
for returns, short-time dependence in the autocorrelation 
function of returns and, on the contrary, long-term de- 
pendence in the volatility signal which shows the typical 
clustering and "declustering" of the real market data. 
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ity dataset of the change rate EURO /USD and from the 
stock price p(t) obtained from the Ising model. In Fig- 
ure [7| we report the PDF P(H) and the cumulative PDF 
P(H > h) of the variable H. From both panels, we can 
note a very similar statistical behaviour and in partic- 
ular, a clear departure from a local Poisson hypothesis 
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